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Abstract-  An  analytical  method  is  presented  to  investigate 
the  stresses  and  strains  in  the  cartilage  layers  of  diarthrodial 
joint  under  rolling  contact.  Each  cartilage  layer  in  contact  is 
assumed  to  be  biphasic,  composed  of  a  linear  elastic  solid  phase 
and  a  Newtonian  viscous  fluid,  whereas  the  subchondral  bone 
is  simplified  as  a  rigid  body.  The  contact  range  of  two  cartilage 
layers  is  discretized,  where  the  surface  tractions  are  piecewise 
given.  A  Galerkin-penalty  method  is  applied  to  form  the  finite 
element  formulation  for  the  cartilage  layers.  The  surface 
tractions,  stresses  and  strains  in  the  cartilage  layers  are  then 
obtained  for  each  time  step  with  a  numerical  procedure  of 
rolling  contact.  Results  show  that  the  interstitial  fluid  plays  a 
fundamental  role  in  the  distributions  of  the  stresses  and  strains 
in  the  cartilage.  The  normal  solid  stress  reaches  its  maximum 
on  the  cartilage-subchondral  bone  interface.  The  coefficient  of 
friction  at  the  contact  surface  has  a  great  effect  on  the 
tangential  traction  while  it  has  little  effect  on  the  normal 
traction.  The  difference  of  rolling  velocity  between  the  two 
cartilage  layers  has  an  increasing  effect  on  the  tangential 
traction  as  the  coefficient  of  friction  increases. 

Keywords-  Rolling  contact.  Articular  cartilage,  Diarthrodial 
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I.  INTRODUCTION 

The  function  of  diarthrodial  joint  is  to  transmit  loads  and 
allow  motion  between  the  bones  of  the  musculoskeletal 
system.  It  is  considered  that  the  mechanical  failure  of 
articular  cartilage  in  diarthrodial  joint,  such  as  keen  joint,  is 
the  direct  cause  of  osteoarthritis.  A  variety  of  theoretical 
models  have  been  developed  for  estimating  the  contact 
pressure  or  the  stress  distribution  in  loaded  articular 
cartilage  [  1  ]~[3] .  However,  it  is  necessary  to  make  further 
research  considering  various  factors,  including  the  friction 
of  two  cartilage  layers  under  rolling  contact,  for  overall 
understanding  of  the  mechanism  of  degeneration. 

In  this  paper  an  analytical  model  of  two  cartilage  layers 
under  rolling  contact,  which  is  based  on  the  KLM  linear 
biphasic  theory,  is  developed.  A  piecewise  discretization 
method  together  with  a  FEM  procedure  is  adopted  to  obtain 
the  mechanical  response  of  articular  cartilage. 

11.  METHODS 

A.  Contact  model 

Figure  1  shows  the  contact  configuration  of  the  two 
cylindrical  biphasic  cartilage  layers  of  radii  R{,  thickness  bv 
and  surface  velocities  V{  (7=1,  2)  with  a  contact  width  2 a  and 
a  normal  approach  d  under  a  normal  load  per  unit  length  P 
along  the  longitudinal  axis  of  the  cylinder.  The  subchondral 
bone  of  each  layer  is  assumed  impermeable  and  rigid  in  this 
analysis.  The  inertial  force  in  each  layer  is  ignored  as  it  is 


sufficiently  small  compared  to  the  contact  pressure.  Each 
cylindrical  layer  can  be  approximated  by  a  flat  layer  bonded 
to  a  rigid  half-plane  if  the  contact  width  is  small  compared 
to  the  cylinders'  radii. 


B.  Governing  equations 

According  to  the  linear  biphasic  theory  [4],  the  solid  and 
fluid  stresses,  cr'  and  af,  are  related  to  the  interstitial  fluid 
pressure  p  and  the  effective  solid  stress  O'  by 

a  =  -<p'pl  +  a"  (la) 

af  —  —tpfp\  (lb) 

where  (p  is  the  volume  fraction  (<ps  +  <pf  =  1),  A  and  p  are  the 
Lame  constant  of  the  solid  phase,  which  together  define  the 
aggregate  modulus  HA  =A  +  2p.  The  effective  stress  G  is 
relate  to  the  solid  strain  e  in  the  form: 


o'  =  Atr(e)l  +  2Jue 

(2) 

The  solid  and  fluid  quasistatic  momentum  equations,  and  the 
continuity  equation  are  given  by 

Vos  +K(\f  -  vs)  =  0 

(3a) 

< 

Cl 

1 

< 

i 

ii 

o 

(3b) 

V((ps\s  +  (pf\f)  =  0 

(3c) 

where  v  is  the  phase  velocity,  and  the  drag  constant  K  is 
related  to  the  fissure  permeability  k  by  K  ={tpf)2l  k. 


C.  Boundary’  conditions 

Define  the  relative  fluid  flux  w  =(ff{y/-  Vs)  and  the  solid 
displacement  u.  Since  the  subchondral  bone  of  each  layer  is 
assumed  impermeable  and  rigid  in  this  study,  the  cartilage- 
bone  boundary  conditions  are  as  follows: 

At  y  =  b ,  (ux)i  =  ( uy)t  =  (wy)i  =  0 ,  -°°  <  x  <oo  (4) 

At  the  contact  interface,  the  boundary  conditions  stipulate 
continuity  of  the  fluid  pressure,  effective  traction  normal 
displacements  and  normal  fluid  flux  within  the  contact 
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region.  Thus: 

At  y  =  0,  - p,  +  (a;,),.  =  ty  (x) ,  (<,),.  =  tx (x)  (5a) 

2 

£(«„),.  =  d  +  ex/R-x2/2R,  H<a  (5b) 

1=1 

LK),  =  o  (5c) 

i= 1 

where  tx(x),  t(x)  are  the  normal  and  tangential  traction  at  the 
contact  interface  respectively,  e  is  an  offset  between  the 
contact  center  and  the  line  which  connects  the  centers  of  the 
two  cylinders  ,  and  R  is  the  equivalent  radius  of  curvature, 
R1R2/(Ri+R2)-  Outside  of  the  contact  region  the  surface 
tractions  reduce  to  zero,  so  that: 

At  y  =  0,  p,  =  0  ,  (&>„),  =  0  ,  «.),  =  0 ,  |x|  >  a  (6) 

Finally,  the  total  load  per  unit  length  P  applied  across  the 
surface  must  be  balanced  by  the  normal  traction  at  the 
contact  interface, 

a 

P=  jty(x)dx  (7) 


In  a  steady-state  of  rolling  contact,  the  surface  velocities 
f  (i=  1,  2)  of  the  layers  1  and  2  in  the  contact  region  can  be 
written  as  follows: 


f=VX  1 +  (£,),],  |x|  -  a  (8) 

The  contact  region  is  divided  into  two  parts:  one  is 
called  the  stick  zone  where  two  contacting  layers  have  the 
same  surface  velocity,  and  another  is  called  the  micro-slip 
zone  where  the  two  layer  surfaces  slid  each  other.  At  the 
stick  zone,  we  have  from  (8): 

l^l  +  fe)^  (9) 

U2  1 +  (£,), 

In  the  micro-slip  zone,  the  Coulomb's  law  is  assumed  to 
satisfy  and  then 

tx  (x)  =  sgn(/2  -  f  )vt  (x)  ( 1 0) 


where  v  is  the  coefficient  of  friction  between  two  layer 
surfaces. 


D.  Solution 

The  contact  region  (-a,  a)  shown  in  Fig.l  is  equally 
divided  into  2 S  elements  and  the  distributions  of  the  normal 
and  tangential  traction,  tx(x),  t,(x),  in  the  contact  region  are 
approximated  by  a  series  of  straight  lines  defined  by  txn)  and 
tyn)  at  the  points  of  intersection,  n,  as  shown  in  Fig. 2.  These 
approximations  are  equivalent  to  that  tx(x)  and  ty(x)  are  each 
represented  by  25-1  overlapping  triangle  elements  with 
equal  base  of  width  2 a/S.  The  tractions  represented  by  the 
nth  triangle  element  in  the  local  coordinate  system  (x*,  y*) 
originated  at  the  point  n  are  given  in  the  form: 

4 (x*)  =  txn) (l  -  S / a\x  *|)  (11a) 

Fv(x*)  =  4n)(l-5/a|x*|)  (lib) 

A  Galerkin-penalty  finite  element  method  [5]  is  then 
employed  to  solve  the  displacements  and  stresses  of  each 
layer  under  the  tractions  shown  in  Fig.  2.  In  this  approach,  a 
weighted  residual  statement  is  constructed  for  a  solid  phase 


and  for  the  fluid  phase,  and  the  penalty  method  is  used  to 
introduce  the  continuity  equation.  The  governing  first  order 


Fig.  2.  Piecewise  discretization  of  contact  region 


coupled  differential  equations  for  the  finite  element 
assembled  system  are  in  the  form: 

Cv  +  Ku=F  (12) 

where  C,  K  and  F  are  the  conductivity  matrix,  stiffness 
matrix  and  force  vector,  respectively.  We  solve  (11)  using 
finite  difference  techniques.  Considering  times  t„,  and  tm+l 
separated  by  a  time  increment  A  t  =  tm+l  -  t„,  and  using  the 
Generalized  Trapezoidal  Family  of  first  order  difference 
rules,  We  obtain: 

[C+  co  Af  K  ]v„,+1  =  F„,+1  -  K  K+A?(l-  co)v„], 

0  <  co  <  1  (13) 

Equation  (12)  is  a  velocity  based  formulation  and  can  be 
solved  recursively  for  vm+1  once  a  set  of  boundary  conditions 
and  initial  conditions  have  been  given.  Then  the  nodal 
displacements  um+1  can  be  computed  by 

u„,+i=u„,+  At[(l-co)v,„+covm+1]  (14) 

We  can  obtain  the  displacements  of  the  contact  surface 
in  each  layer  under  the  «th  traction  element  in  which  txn)  and 
tyn)  («=-5+l ,  •••,  5-1)  are  put  to  unit  in  (11a)  and  (lib)  by 
using  (13)  and  (14).  According  to  (9),  (10)  and  the  above 
results,  we  can  solve  the  unknowns  txn)  and  tyn)  (n=-S+ 1,  •••, 
5-1)  when  the  contact  half- width  a  or  normal  approach  d, 
surface  velocity  1A  (i=  1,  2)  and  so  on  are  given.  Then 
solving  (13)  and  (14)  again,  we  can  calculate  the 
displacements  u  and  phase  velocities  v  due  to  the  whole 
surface  tractions  shown  in  Fig. 2.  The  interstitial  fluid 
pressure  in  the  cartilage  layers  can  be  obtained  by 

p  =  -j5  V(0V  +  0V)  (15) 

where  (5  is  the  penalty  parameter  which  is  selected  large 
enough  to  enforce  the  continuity  constraint  condition,  (3c), 
but  not  so  large  that  the  governing  matrix  equations  become 
ill  conditioned.  Finally,  we  obtain  <r5,  Of  and  &’  according  to 
( la),  (lb)  and  (2). 

111.  RESULTS  AND  CONCLUSIONS 

In  the  present  results,  a=2mm,  bl=b2= 2mm,  K,=8mm, 
R2=6mm,  A=0.05MPa,  R=0.25MPa,  <p- 0.2,  (ft= 0.8,  /~4.0  x 
10‘15m4/Ns. 

Figure  3  shows  the  contact  tractions  for  two  values  of 
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the  coefficient  of  friction  v  (0.001  and  0.3).  It  is  seen  that  v 
has  a  great  effect  on  the  tangential  traction  while  it  has  little 
effect  on  the  normal  traction. 

Figure  4  shows  the  distribution  of  the  normal  stress  for 
y/b,= 0.05,  0.55  and  0.95  in  layer  1.  It  is  seen  that  the  normal 
stress  reaches  its  maximum  on  the  cartilage-subchondral 
bone  interface. 


xta 

Fig.  4.  Effective  solid  Stress  O'jj,  fory6,=0.05,  0.55  and  0.95 
in  layer  1(F,/F2=0.95) 


xta 

Fig.  5.  Interstitial  fluid  pressure p  lot'  t-.  /r ,  0.05 .  0.55  and  0.95 
in  layer  1(F,/K2=0.95) 


Figure  5  shows  the  interstitial  fluid  pressure  p  for 
ylbx= 0.05,  0.55  and  0.95  in  layer  1.  It  is  observed  that  the 
interstitial  fluid  pressurization  could  contribute  more  than 
80%  of  the  total  applied  load  at  the  contact  surface. 

Figure  6  displays  the  tangential  traction  tx(x)  for 
Vx/V 2=0.95  and  0.90  in  layer  1.  It  is  observed  that  the 
magnitude  of  the  tangential  traction  increases  as  the 
difference  of  rolling  velocity  between  the  two  cartilage 
layers  enlarges.  Besides,  the  effect  of  the  difference  on  the 
tangential  traction  becomes  larger  as  the  coefficient  of 
friction  increases. 


xla 

Fig.  6.  Tangential  traction  tx(x)  for  Fi/F2=0.95  and  0.90  in  layer  1 
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